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Abstract 

The Pierre Auger Collaboration has reported evidence for anisotropy in the distribution 
of arrival directions of the cosmic rays with energies E> Eth = 5.5 X 10" eV. These show a 
correlation with the distribution of nearby extragalactic objects, including an apparent excess 
around the direction of Centaurus A. If the particles responsible for these excesses bX, E > Eth 
are heavy nuclei with charge Z, the proton component of the sources should lead to excesses in 
the same regions at energies E/Z. We here report the lack of anisotropics in these directions 
at energies above EthjZ (for illustrative values of Z = 6, 13, 26). If the anisotropies above 
Eth are due to nuclei with charge Z, and under reasonable assumptions about the acceleration 
process, these observations imply stringent constraints on the allowed proton fraction at the 
lower energies. 
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1 Introduction 



Anisotropy and composition, together with the study of the features in the energy spectrum, are 
the fundamental tools available to decipher the origin and nature of the ultra-high energy cosmic 
rays (UHECRs). The suppression of the flux observed above 40 EeV suggests that the energy 
of the UHECRs is attenuated by interactions with the cosmic microwave background and infrared 
photons on their journey from their extragalactic sources, either by photopion interactions in the 
case of protons or by photodisintegration in the case of nuclei i4| . This would imply that at 
the highest energies cosmic rays can only arrive from nearby sources, within the so-called GZK 
horizon (which is e.g. ^ 200 Mpc for protons above 60 EeV O 1^). This is supported by the 
correlation reported by the Pierre Auger Collaboration O IHl |^ between the arrival directions 
of cosmic rays with energies above 55 EeV and the distribution of nearby extragalactic objects. 
The correlation with nearby active galactic nuclei (AGN) in the Veron-Cetty and Veron (VCV) 
catalog [10 was originally found with data collected up to May 2006, and was most significant 
for the AGN within 75 Mpc and for angular separations between the AGN and cosmic ray arrival 
directions smaller than 3.1°. A test with subsequent data rejected the null hypothesis of isotropy 
with 99% confidence [71 [8]. A more recent analysis [9 has found that the fraction of events above 
55 EeV correlating with these AGN is (38lg)%, smaller than obtained initially but still well above 
the isotropic expectation of 21%. Note that these AGN may well be acting just as tracers of the 
actual UHECR sources, and indeed it is interesting that alternative studies with other populations 
(X-ray AGN from the SWIFT catalog or galaxies from the 2MASS catalog) also indicate some 
degree of correlation within a few degrees with those objects l9l (see also [HI [12j [131 (HI US])- 
The final identification of the UHECR sources will require much additional data. 

Another interesting cosmic ray excess was found in the direction towards Cen A, at equatorial 
coordinates {a, 6) = (201.4°, —43.0°). Already in [71 [8] it was pointed out that two out of the 
27 highest energy events observed before August 2007 by the Pierre Auger Observatory arrived 
within less than 3° of Cen A, with several more events lying in the vicinity of its radio lobes. 
More recently, with data up to the end of 2009 9 and considering the events above 55 EeV, the 
most significant excess around Cen A was identified for an 18° window, in which 13 events were 
observed while only 3.2 were expected. Whether this excess, if confirmed with further data, is due 
to Cen A, which is one of the nearest AGN (being at less than 4 Mpc distance), or due to one 
or several sources farther away, e.g. in the Centaurus cluster lying in a similar direction but at 
^ 45 Mpc, is something that remains to be determined. It should be mentioned that the HiRes 
air shower experiment has not found indications of an excess correlation with nearby AGN |17] . 
although the associated statistics are smaller and there are systematic differences in the energy 
calibrations between the two experiments. Also, contrary to Auger, the HiRes experiment looks 
to the northern hemisphere, and in particular this makes it blind to the Cen A region of the sky. 

The Pierre Auger Observatory has recently measured the average depth of the maximum of 
shower development X^ax and its fiuctuations |18j . The logarithmic slope of the average shower 
maximum vs. energy becomes smaller above ~ 2 EeV, indicating a change in the shower properties. 
Also the fiuctuations in X^ax become suppressed above this energy. An inference of the chemical 
composition of the primary cosmic rays can be done via comparison with Monte Carlo simulations 
of air showers. If these models are taken at face value, they indicate a gradual increase in the 
average mass as a function of energj0. Alternatively, this behavior could be ascribed to changes 
in the hadronic interactions (cross sections, inelasticities or multiplicities) not considered in the 
available models. We note that the models make extrapolations to energies well beyond those 
tested at accelerators. One should also keep in mind that, due to the limited statistics of the 

^ We note that the HiRes experiment measures a depth of shower maximum consistent with proton-only Monte 
Carlo air shower simulations all the way from 1 EeV up to ~ 40 EeV |19| . 
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events observed with fluorescence telescopes, there is no measurement of the mean Xmax and the 
corresponding fluctuations available for > 55 EeV. It is clear that performing alternative studies 
to try to improve our understanding of the UHECR composition is important. 

In this work we perform searches for anisotropics in the same directions where excesses were 
observed above Eth = 55 EeV, but using lower energy thresholds (we consider the illustrative 
values for the threshold Eth/Z, with Z = 6, 13, 26). 

We first focus on the analysis of the region around Cen A, for which the most significant excess 
was found above 55 EeV for an angular window of 18° radius. We note that the location of the 
excess, the size of the angular window and the selected energy threshold are a posteriori, therefore 
new independent data would be required to assess the significance of the excess at high energies. 
However, we are already able to report the results of our search for anisotropics in the same region 
for the lower energy thresholds considered. 

We also perform a similar search but looking for possible excesses in windows of 3.1° around 
the VCV AGN within 75 Mpc, considering only the data after May 2006 so as to exclude those 
used to fix these parameters. We note that even above 55 EeV the cosmic ray deflections in the 
galactic magnetic field are likely larger than a few degrees (especially if cosmic rays happen to be 
heavy nuclei). However, the VCV correlation does not imply that the objects in this catalog are 
the sources, nor that the typical deflections are smaller than the optimal correlation angle. Active 
galaxies in the VCV catalog trace the nearby large scale matter distribution, and that includes all 
types of candidate astrophysical sources, not only AGN and their subclasses. Defiections of cosmic 
ray trajectories could be larger, and still manifest an anisotropy through a correlation of a fraction 
of them within a few degrees of the VCV objects. 

We then explore the possibility that the anisotropics at the highest energies might be due to 
heavy nuclei. Using our observations and following an idea proposed by Lemoine and Waxman 
[20] , which exploits the fact that a high energy anisotropy due to nuclei of charge Z should lead to 
an anisotropy in the same region of the sky at energies Z times smaller due to the protons from 
the same sources, we are then able to constrain the allowed proton fraction at the source under 
different assumptions on the value of the nuclear charges responsible for the high energy excess. 

2 The Observatory and the dataset: 

The Pierre Auger Observatory is located near the town of Malargiie, Argentina, at a latitude of 
35.25° S. It is a hybrid detector, consisting of 24 fluorescence telescopes and a surface array of 
1600 water Cherenkov detectors covering ^ 3000 km^ (see (21] [22] for further details). 

The data considered in the present work consists of the cosmic ray events with zenith angles 
9 < 60° detected by the surface array (which has an almost 100% duty cycle and hence collected 
the largest data set) since 1 January 2004 up to 31 December 2009. The array has been growing 
in size until the completion of the baseline design in mid 2008. In order to have an accurate 
estimate of the exposure and hence of the expected background in the different regions of the sky 
we have removed periods in which the data acquisition was unstable (the resulting livetime being 
87% [21]) and applied a quality cut that requires that for any event the six detectors surrounding 
the detector having the largest signal be active at the time the event is recorded. Keeping track 
of the number of active detector configurations able to trigger such events at any time allows us 
to take into account the detector growth and dead times in the evaluation of the exposure. The 
isotropic expectation in an angular window Afl can be obtained as iV^so = x Ntot, where x is the 
fraction of the exposure within the solid angle Ail and Ntot is the total number of events. 

The trigger efficency is 100% for i? > 3 EeV, but at lower energies (we consider here events 
down to £^ = 55 EeV/26 ~ 2.1 EeV) the trigger efficiency becomes smaller than unity and is 
zenith angle dependent. Hence, to obtain the isotropic expectations for the lower energy threshold 
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considered we use a fit to the zenith angle distribution of the events, rather than the ideal exposure 
expectation dN ex. sin 9 cos 9dd. We note that the detection efhciency below 3 EeV may also depend 
on the composition of the cosmic rays, being actually smaller for lighter nuclei. This could slightly 
affect the predictions for the expected localized proton excesses for the lowest energy threshold 
considered, E > 2.1 EeV. We estimate that in this case the predictions are affected by no more 
than 2% by the possible differences in exposure (using the values in ref. [21]), and hence these 
effects can be safely neglected. 

3 Results 

3.1 The Centaurus A excess 

We first consider the excess observed in the Cen A region for energies above the threshold Eth — 
55 EeV. The cumulative number of events as a function of the angular distance from the direction 
of Cen A is plotted in Fig. [T] (to make the plot more readable we display the difference with respect 
to the average isotropic expectations) . 

10 I 1 1 1 1 1 1 




Angular distance to Cen A (degrees) 

Figure 1: Cumulative number of events with E > 55 EeV (subtracting the average isotropic 
expectations) as a function of angular distance from the direction of Cen A. The bands correspond 
to the 68%, 95% and 99.7% dispersion expected for an isotropic flux. 

In Fig. [2] we plot the cumulative number of events, subtracting the isotropic expectations, 
as a function of the angular distance from the direction of Cen A for lower energy thresholds, 
considering energies above Eth/Z in the cases Z = 6, 13 and 26. The observed distributions are 
consistent with the isotropic expectations (shaded regions), showing no significant excesses in any 
of the angular windows considered. 

As reported in [9J , the most significant excess for a top-hat window around Cen A was obtained 
for a radius 7 = 18° and we will hence focus on this region. For this energy range, the total number 
of event^ is Ntot — 60, with Nobs = 10 of these being in an 18° angular window around Cen A. 
If we adopt the expression for the ideal exposure of the detector, the fraction of isotropic sky in 

^ Different from ref. f9l, where 13 out of 69 events were reported to correlate within 18° of Cen A, the stricter 
event selection applied in this work in order to get an accurate estimate of the exposure at low energies yields 10 
correlations out of 60 events, well within the statistical uncertainties of the previous result. 
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this 18° region is x ~ 0.0466. Normalizing to the counts outside the source region, the expected 
background in this region is N^kg — {Ntot — Nobs)x/{l — x) — 2.44 counts. 





Angular distance to Cen A (degrees) Angular distance to Cen A (degrees) 



Figure 2: Similar to Fig. [TJ events with E > 55 EeV/Z for Z = 6 (top), 13 (bottom left) and 26 
(bottom right). 



In table[T]we report the observed number of events with E > 55 EeV/Z (total and in an angular 
window of 18° around Cen A), as well as the expected isotropic background. No significant excess 
is found for any of the lower energy thresholds considered. 



z 


E,mn [EeV] 


Ntot 


Nobs 


Nbkg 


6 


9.2 


4455 


219 


207 ± 14 


13 


4.2 


16640 


797 


774 ± 28 


26 


2.1 


63600 


2887 


2920 ± 54 



Table 1: Total number of events. Not, and those observed in an angular window of 18° around 
Cen A, Nobs, as well as the expected background Nbkg- Results are given for different energy 
thresholds, corresponding to Emin — Etu/Z for the indicated values of Z and Eth = 55 EeV. 



3.2 The VCV AGN 

We now search for possible overdensities of cosmic rays with arrival directions within 3.1° of 
objects with redshift z < 0.018 (~ 75 Mpc) in the VCV catalog. We use for this study only data 
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collected after May 2006, subsequent to data used to specify the parameters that optimized the 
VCV correlation in that period. 

In this case, one has that for E > Eth — 55 EeV there are Ntot = 49 events, of which Nobs = 20 
are within 3.1° of the nearby AGN. On the other hand, the probability that isotropic cosmic rays 
correlate by chance with those objects is cc ~ 0.212 and hence N^kg = {Ntot — Nobs)x/ (l — x) = 7.88. 

In table Owe show the observed number of events with _E > 55 EeV/Z (total and those within 
3.1° of an object with z < 0.018 in the VCV catalog), as well as the expected background. It is 
apparent that no significant excess is found for any of the lower energy thresholds considered. 



z 


S™„ [EeV] 


Ntot 


Nobs 


Nbkg 


6 


9.2 


3626 


763 


770 ± 28 


13 


4.2 


13482 


2852 


2860 ± 54 


26 


2.1 


51641 


10881 


10966 ± 105 



Table 2: Total number of events, Ntot, and those observed within 3.1° from objects with z < 0.018 
in the VCV catalog, Nobs, as well as the expected isotropic background N^kg- Results are given 
for different energy thresholds, corresponding to Emin = Eth/Z for the indicated values of Z and 
Eth = 55 EeV. 



4 Constraints on the source composition 

As a by-product of the observations described above, and under reasonable assumptions on the 
cosmic ray acceleration and propagation, it is possible to set some constraints on the composition 
of the cosmic rays responsible for localized overdensities observed above Eth- In order to do this, 
we elaborate on an idea proposed by Lemoine and Waxman |20j . who related the high energy 
excess, under the assumption that it is due to heavy nuclei of charge Z, with the expected excess 
at energies above Eth/Z due to the protons from the same sources. Note that, in the absence of 
energy losses and scattering effects, protons with energies E/Z would follow the same trajectories 
as nuclei of charge Z and energy E coming from the same source, and hence they should arrive 
within the same angular windows. Moreover, even if at lower energies the isotropic background 
can be enhanced by the contribution from sources beyond the GZK horizon, the gain in statistics 
obtained can make the search sensitive to relatively small low energy anisotropics. 

The main underlying hypothesis is that the cosmic ray acceleration depends just on the particle 
rigidities, i.e. on E/Z. It is therefore natural to assume that at the sources the spectra of the 
different charge components scale as 

^^kzm/z), (1) 

with kz being constant factors. The function $ may display a high energy cutoff resulting from 
the maximum rigidities attainable by the acceleration process. If, in this scenario, the maximum 
proton energies were below Eth, the higher energy cosmic rays from the source could be dominated 
by a heavy component. 

If A^(> E) is the number of events with energies above the threshold E which come within a 
certain solid angle around a source and if the acceleration process at the source depends only on 
rigidity, then the number of nuclei of charge Z above Eth and those of protons above Eth/Z are 
related by 

Afp(> Eth/Z) = -^Nz{> Eth). (2) 
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This relation does not take energy losses into account (included as a parameter a in ref. pU]). 
Ignoring them leads to more conservative bounds on the ratio kp/kz, because energy losses are 
larger for nuclei of charge Z and energy E than for protons of energy E/Z. Moreover, the nucleons 
emitted in the photodisintegration processes can also add to the expected proton anisotropics at 
low energies. 

The number of events produced by the source(s) responsible for the localized excess observed 
can be estimated as TV = -/Vof,s — -^bfcg in terms of the number of events observed in the window 
considered and the expected background, which are displayed in the tables. Taking into account 
the Poisson fluctuations in the low and high energy signals, as well as in the background estimates, 
we obtain 95% CL upperbounds on the quantity Rz = N{> Eth/Z)/N{> Eth) using the profile 
likelihood method (see e.g. [13]). In the case of Cen A these bounds are Rz < 12.9, 17.3 and 9.1 
for Z = 26, 13 and 6 respectively, while for the case of VCV the bounds are Rz < 14.7, 12.4 and 
6.0. We note that considering a 99% confidence level, the bounds become Rz < 23.8, 31.1 and 
16.3 for Cen A, and Rz < 28.9, 23.7 and 11.4 for VCV (for Z = 26, 13 and 6 respectively), being 
then typically a factor of two weaker. 

If the excess at high energies is indeed dominated by the heavier nuclear component of charge 
Z, i.e. iV(> Eth) ^ Nz{> Eth), we obtain that Rz > Np{> Eth/Z)/Nz{> Eth) + l = kp/{Zkz) + l 
(where we used that Nz{> Eth/Z)/Nz{> Eth) > 1)- In this way, conservative bounds kp/kz < 
Z{Rz — 1) can be obtained. 

One may translate these limits on the relative spectrum normalizations into bounds on the 
actual low energy abundance ratios between the proton and heavy elements at the source. In 
particular, in the case that one assumes that below a certain rigidity the spectrum has a power 
law behavior, i.e. $ cx {E/Z)~'' for E/Z < Ei, as expected in scenarios of diffusive shock ac- 
celeration, at energies below Ei all the relative abundances of the different elements present will 
be independent of the energy. In this case, one can relate the low energy relative fractions ft of 
the different elements at the source with the normalization factors kt in eq. ([T]). Comparing the 
differential spectra for protons and for the charge Z at energies below Ei one gets 

kp _ fp „s 

Note that we are not making any assumption about the spectral shape above the threshold energies. 
Also, for energies above Ei the values of fp/fz will depend on the spectral shape details. 

The resulting bounds for the low energy relative abundances are displayed in Figs. [3] and S] 
as a function of the low energy spectral index, for values 1.5 < s < 2.5. The regions above the 
respective lines are excluded at 95% CL. 

In the case in which energy losses can be neglected, such as if the source of the excess events 
is the nearby Cen A galaxy, it is appropriate to consider energy bins and relate, through an 
expression analogous to that in eq. ^ , the events in the bin [Eth , "^Eth] to those at Z times lower 
energies (where we adopted for definiteness a bin width corresponding to a factor two in energy 
The ratio between the events observed and those expected for the 18° window around Cen A are 
Nobs/Nbkg = 152/153.5, 543/533.4 and 2090/2147.9 for Z = 6, 13 and 26 respectively. This leads 
to bounds on the proton fractions similar to those in Fig. |3]but about a factor two stronger. 

An important point is that the statistical significance of the constraints in Fig. [3] is a posteriori, 
since the identification of the region around Cen A, its angular size and the energy threshold 
were tuned to maximize this excess. Therefore it would be necessary to look to this same region 
using the same energy threshold with an independent dataset of comparable size so as to obtain an 



•^If energy losses were relevant, the observed energies of the high energy events might correspond to a wider 
span of energies at the source, and the corresponding low energy protons may then span a range of energies wider 
than [Eii^/Z,2Eti-^/Z], making this analysis no longer valid, while that based on the integral energy bins above a 
threshold would still provide conservative bounds. 
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Figure 3: Upper bounds at 95%CL on the allowed proton to heavy fractions in the source as a 
function of the assumed low energy spectral index s. The different lines are for charges Z = 6, 13 
and 26, as indicated. 



unbiased estimate of the strength of the source (or sources) producing the excess. We note however 
that varying the energy threshold to 50 or 60 EeV leads to qualitatively similar results. Also the 
angular size adopted for the window is not very crucial. For instance, if we consider a 10° window 
instead of the 18° one, the main effect on the bounds comes from the modification of the expected 
background in the low energy bin. This would relax the bounds on fp/ fpe by a factor of about 
two in this case. It is interesting to point out that the sensitivity of the a posteriori bounds from 
Cen A turns out to be comparable to that achieved with the analysis of the VCV correlations. 



5 Discussion and conclusions 

We have searched for overdensities at energies Eth/Z in the regions where anisotropics were re- 
ported previously above Eth = 55 EeV, i.e. both in the direction towards Cen A and in 3.1° 
windows around nearby AGN from the VCV catalog. Considering representative values of Z = 6, 
13 and 26, where Z is the assumed charge of the cosmic rays responsible for the high energy 
anisotropics, we have found no indications of overdensities in any of the lower energy bins. In 
scenarios where the acceleration process is only dependent on rigidity, the absence of significant 
anisotropics at energies E > Eth/Z implies that an upper bound can be set on the low-energy 
relative proton abundance at the sources. This bound is given by fp/fz < (0.5 to 2)Z^^'' at the 
95% CL level, depending on the adopted value of Z (see lines in Figs. [3] and 2]). Note that the 
constraints become weaker if the source spectrum is very hard (s ~ 1.5). Given the comparable 
bounds obtained for different values of Z, similar limits will result in the case in which the high 
energy anisotropy is dominated by nuclei belonging to a given mass group with similar values of 
Z. 

On the other hand, estimates of the expected low energy relative abundances point towards 
values above these bounds. For instance, the ATIC-2 experiment [21] has measured that at 100 TeV 
(the highest energies for which detailed composition measurements are available) one has /p ~ 
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95% CL upper bounds from VCV 



N 

a. 




0.1 I 1 ' ' 1 ' ^ 

1.6 1.8 2 2.2 2.4 

S 

Figure 4: Upper bounds at 95% CL on the allowed proton to heavy fractions in the source as a 
function of the assumed low energy spectral index s implied by the absence of an excess of arrival 
directions of cosmic rays with energies above 55 EeV/Z within 3.1° of objects with z < 0.018 in 
the VCV catalog. The different lines are for charges Z — 6, 13 and 26, as indicated. 



Jhc — "^fcNO — — '^fz>i7 — ^fpe- Morcover, for these low energies (for which cosmic 

rays are believed to be of galactic origin) one would expect that the measured relative fraction 
of protons versus heavy nuclei for a given particle energy is actually smaller than the original 
fraction at the sources, due to the longer confinement time in the Galaxy of the heavier species. 
For instance, the measured p to Fe fraction would be 26^/'^ ~ 3 times smaller than the value at the 
source if the turbulent component of the galactic magnetic field has a Kolmogorov spectrum, so 
that the diffusion coefficient scales as Z? cx [E/ZY/'^ . One has to keep in mind that these estimates 
based on lower energy galactic cosmic ray sources do not necessarily apply to the extragalactic 
sources which are most likely responsible for the highest energy events, but one may consider that 
they provide a useful indication of the plausible expected values. 

Hence, we conclude that a heavy composition for the excesses observed at high energies appears 
to be in conflict with rigidity-dependent acceleration scenarios having at low energies a proton 
component more abundant than heavier species, as quantified in Fig. 3] How these conclusions are 
modified in the presence of strong structured magnetic fields and taking into account the relevant 
energy losses remains to be seen. We note that the present analysis based on the lack of anisotropics 
at lower energies provides information on the cosmic ray composition which is independent of Xmax 
measurements, but depends instead on assumptions related to source properties. 
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